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Abstract 

We report measurements of radiative B decays with Krj^ final states, using a data sam- 
ple of 140 fb"^ recorded at the T(45) resonance with the Belle detector at the KEKB e'^e 
asymmetric energy collider. We observe B Kr}^ for the first time with a branching frac- 
tion of (6.9li:^(stat)ti:g(syst)) x lO'^ for Mkt^ < 2.4 GeV/c^ We also set an upper limit on 
B K|(1780j7. 

PACS numbers: 13.20.He, 14.40.Nd 
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Radiative B decays, which proceed mainly through the 6 — S7 process have played 
an important role in the search for physics beyond the Standard Model (SM). Although 
the inclusive branching fraction has been measured to be (3.3 ± 0.4) x 10~"^ j2|, we know 
little about its constituents. So far, measured exclusive final states such as ii'* (892)7 0,01, 
/r* (1430)7 B li, K-KTc-f \E\ and B K(f)'y § only explain one third of the inclusive 
rate. Detailed knowledge of exclusive final states reduces the theoretical uncertainty in the 
measurement of the inclusive branching fraction using the pseudo-reconstruction technique, 
as well as in the measurement of B Xsi'^i~ . In this analysis, the decay mode B — >■ Krj'j 
is studied for the first time. In addition to improving the understanding of 6 — *• 57 final 
states, B^ K^Vl can be used to study time-dependent CP asymmetry 0], which is 
sensitive to physics beyond the SM. The mode B —>■ Krj'-f can also be used to search for 
B is:* (1780)7 via 7^1(1780) ^ Kr] decay. 

The analysis is based on 140 fb~^ of data taken at the T(4S') resonance (on- resonance) and 
15 fb"^ at an energy 60 MeV below the resonance (off-resonance), which were recorded by 
the Belle detector [y] at the KEKB asymmetric e~^e~ collider (3.5 GeV on 8 GeV) The 
on-resonance data corresponds to 152 million BB events. The Belle detector has a three- 
layer silicon vertex detector, a 50-layer central drift chamber (CDC), an array of aerogel 
Cherenkov counters (ACC), time-of-flight scintillation counters (TOF) and an electromag- 
netic calorimeter of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An instrumented iron flux-return for Kl/ ^ detection 
is located outside the coil. 

We reconstruct B'^ — > K'^rj'y and B^ K^rj'y via 77 ^ 77 and f] tt'^ti^ti^. All 
the charged tracks used in the reconstruction (except charged pions from Kg) are required 
to have an impact parameter within ±5 cm of the interaction point along the positron 
beam axis and within 0.5 cm in the transverse plane. Primary charged kaons are also 
required to have a momentum in the e'^e~ center-of-mass (CM) frame that is greater than 
100 MeV/c. In order to identify kaon and pion candidates, we use a likelihood ratio based on 
the light yield in the ACC, the TOF information and the specific ionization measurements 
in the CDC. For the selection applied on the likelihood ratio, we obtain an efficiency (pion 
misidentification probability) of 90% (9%) for charged kaon candidates, and an efficiency 
(kaon misidentification probability) of 98% (9%) for charged pion candidates. 

Kg candidates are formed from tt+tt" combinations whose invariant mass is within 
8 MeV/c^ of the nominal K^ mass. The two pions are required to have a common vertex 
that is displaced from the interaction point. The K^ momentum direction is also required 
to be consistent with the K^ flight direction. Neutral pion candidates are formed from pairs 
of photons that have an invariant mass within 16 MeV/c^ of the nominal 7r° mass and an 
energy greater than 100 MeV in the CM frame. Each photon is required to have an energy 
greater than 50 MeV. A mass constrained fit is then performed to obtain the 7r° momentum. 

1] candidates are reconstructed via — 77 or 77 — >■ 7r+7r~7r°. For t] 77, we require 
that the invariant mass of the two photons satisfy 0.515 GeV/c^ < M^^ < 0.570 GeV/c^ 
and that each photon have an energy greater than 50 MeV. We also require that the r) 
decay helicity angle ^hei satisfies | cos^heil < 0.9. A mass constrained fit is then performed to 
obtain the r] momentum. For rj —>■ tt+tt^tt^, we apply a selection on the Mtt+tt-ttO invariant 
mass, 0.532 GeV/c^ < M^+^-^o < 0.562 GeV/c^. 

We combine a charged or neutral kaon with an rj to form a Krj system with invariant 
mass less than 2.4 GeV/c^. We then reconstruct B meson candidates from the Kt] system 
and the highest energy photon with a CM energy between 1.8 GeV and 3.4 GeV within 
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the acceptance of the barrel ECL (33° < 6^ < 128°, where 9^ is the polar angle of the 
photon in the laboratory frame). The photon candidate is required to be consistent with an 
isolated electromagnetic shower, i.e. 95% of its energy should be concentrated in an array 
of 3 X 3 crystals and no charged tracks should be associated with it. In order to reduce 
the background from decays of 7r° and rj mesons, we combine the photon candidate with all 
other photons in the event and reject the event if the invariant mass of any pair is within 
18 MeV/c^ (32 MeV/c^) of the nominal vr^ {rf) mass (this condition is referred to as the 7r°/?7 
veto). 

We use two independent kinematic variables for the B reconstruction: the beam-energy 

constrained mass Mbc = ^(^Lam/c^/ - (1%, + P^l/c)' and = E*^^ + E; - E*^^^, 
where -Ebeam is the beam energy, and p*, E*, p^^, E^^ are the momenta and energies of the 
photon and the Krj system, respectively, calculated in the CM frame. In the Mbc calculation, 
the photon momentum is rescaled so that \p* \ = (-Ebeam ~ ^Kt]) /c is satisfied. We require 
Mbc > 5.2 GeV/c^ and -150 MeV < < 80 MeV. We define the signal region to be 
Mbc > 5.27 GeV/c^. In the case that multiple candidates are found in the same event, we 
take the candidate that has the rj mass closest to the nominal mass (and smallest |A£'|) 
after applying the background suppression described later. 

The largest source of background originates from continuum e"'"e~ —>■ qq {q = u, d, s, c) 
production including contributions from initial state radiation {e~^e~ —>■ qfl'j). In order to 
suppress this background, we use the likelihood ratio (LR) described in Ref . |^] , which utilizes 
the information from a Fisher discriminant formed from six modified Fox- Wolfram 
moments 2| and the cosine of the angle between the B meson flight direction and the 
beam axis. The LR selection retains 44% of the signal, rejecting 98% of the continuum 
background. 

In order to extract the signal yield, we perform a binned likelihood fit to the Mbc dis- 
tribution. The Mbc distribution of the signal component is modeled by a Crystal Ball line 
shape [l^, where the parameters are determined from the signal MC and calibrated by 
B Dtt decays, as described below. The M^c distribution of the continuum background 
is modeled by an ARGUS function whose shape is determined from the off-resonance 
data. Here, the LR selection is not applied to the off-resonance data in order to compensate 
for the limited amount of data in that sample. The possible bias due to this is taken as 
systematic error to the fitted yield. Background from B decays is divided into two com- 
ponents, which we refer to as BB background and rare B background in this paper. The 
former comprises B decays through b ^ c transitions including color-suppressed B decays 
such as D^TT^, and the latter covers charmless B decays through b ^ s and b ^ u 

transitions. Each of them is modeled by another ARGUS function. The shape of the dis- 
tributions is determined by a corresponding Monte Carlo (MC) sample. In order to study 
the contamination from other 6 — >■ 57 decays, we examine a. B fC* (892)7 MC sample and 
an inclusive b ^ s'y MC sample that is modeled as an equal mixture of sd and su quark 
pairs and is hadronized using JETSET where the Mx^ spectrum is fitted to the model 
of Kagan and Neubert ^16] . We find that the feed-down from other b ^ s'j decays is not 
negligible, and its Mbc distribution is also modeled by an ARGUS function. 

FiguresH(a)-(c) show the Mbc distributions for B~^ —>■ K'^rj'-f, B^ — »• Kgrjj and combined 
B —>■ Krj'-f. The distributions are fitted with signal, continuum, BB, rare B background 
and the b —>■ sj feed-down components. In the fit, the normalization of BB, rare B and 
b ^ sj are fixed according to the luminosity and b ^ s'j branching fraction, while the 
normalization of the continuum component is allowed to float. We flnd the signal yields to 
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FIG. 1: Mbc distributions for (a) (b) B° ^5^7, (c) B Krj-f (combined 

B+ K+rj-f and S° ^ K'^Vl), and (d) S ^ Krij with the mass range for B /r|(1780)7. Fit 
results are overlaid. 



be 35.5 t9°4^, 9.7l;4;3 ^-^d 45.0iJo;4 with statistical significances of 4.3(T, 2.5(T and 5. Oct, for 
the charged, neutral and combined modes, respectively. Here, the significance is defined as 

^—2 ln(£(0)/£max); whcre Cma.x is the maximum of the Ukehhood and C{0) is the hkeUhood 
for zero signal yield. 

The Krj invariant mass distribution for events inside the signal region is shown in Fig. |21 
Here, the background distributions are obtained from the corresponding MC samples, and 
are normalized using the fit result. We find that signal excess is concentrated between 
1.3 GeV/c^ and 1.9 GeV/c^. Therefore, our selection M^ri < 2.4 GeV/c^ is expected to 
include most of the B —>■ Krj'-f signal. We do not see any clear resonant structure in the 
Mxrj distribution. 

The systematic error on the signal yield due to the fitting procedure is estimated by 
varying the value of each fixed parameter by ilex and extracting the new signal yield for 
each case. The difference between the background shape for the continuum MC with and 
without the LR selection is taken as an additional error to the continuum background shape. 
We set the normalization of either the BB or rare B backgrounds to zero and to twice its 
nominal value to account for its uncertainty. The changes of the yields for each procedure 
are added in quadrature, and are regarded as the systematic error on the signal yield. 

The signal reconstruction efficiency is estimated using the MC simulation and is corrected 
for discrepancies between data and MC using control samples. We find that the efficiency 
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is almost independent of the Kt] invariant mass. Table H] shows the signal efficiencies and 
the branching fractions for each B — > Krj'-f mode. Here, we assume an equal production 
rate for B^B^ and B^B~ . The error includes the following systematic uncertainties: photon 
detection (2.8%), tracking (1% track), kaon identification (0.8%), pion identification (0.5% 
per pion), detection (4.5%), 7r° detection (1.5%), rj detection in r/ 77 mode (2.0%), 
7r°/r7 veto and LR (3.2% and 9.0% for charged and neutral modes, respectively), possible 
Krj mass dependence of the efficiency (2.1% and 8.4% for charged and neutral modes, re- 
spectively), and uncertainty in the 77 branching fraction (0.7% for 77 — > 77 and 1.8% for 
7] — s> tt+tt^tt''). The systematic errors from the -k^ /r] veto and LR requirement are estimated 
using the B^ — > D^{-^ K~ti^ti^)tx^ and B^ D^{-^ K^ti^ 71^)11^ as a control sample, 
treating the primary pion as a high energy photon. This sample is also used to obtain the 
Mbc shape of the signal component. 

We also search for the decay B — > fr|(1780)7 by applying the additional requirement 
1.60 GeV/c2 < MKr^ < 1.95 GeV/c^. The ffis to the Mbc distributions yield 10.5i||lli 
4.2 Iq' I In 'o and 15.0l5;5l|;5 events for charged, neutral and combined modes, respec- 
tively 17|. The Mbc distribution and ffi result for the combined mode is shown in Fig. ^(d). 



However, we provide only upper limits due to our inability to distinguish B —>■ 7^^^(1780)7 
from non-resonant decays. The 90% confidence level upper limit is calculated from the 
relation C{n)dn = 0.9 £{n)dn, where C{n) is the maximum likelihood in the Mbc 
fit with the signal yield fixed at n. In order to include the systematic errors from the fit- 
ting procedure in the upper limit for the yield, the positive systematic error is added to 
A^. The obtained yield upper limits, efficiencies and branching fractions are listed in Ta- 
ble H] Here, the error for the efficiency also includes the uncertainty in the i^'3(1780) —>■ Krj 
branching fraction ((30 ± 13)%). The measurement improves the limits set by the ARGUS 
collaboration [Ts^ . 

In conclusion, we observe the decay mode B Krj'y with a branching fraction of 
(6.9iJ;6(stat)lj:o(syst)) x 10"^ for M^r, < 2.4 GeV/c^. We also set upper limits on 
B /Tg (1780)7. Although the signal yield of 5° KgTj'y is small, in future this mode can 
be used to study time-dependent CP asymmetry in radiative B decays and to search for 
new physics. 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
Cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the National Institute of Informatics for valuable computing and Super-SINET network 
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TABLE I: Measured signal yields, efficiencies, branching fractions (B), and statistical significances 
(snf.) Efficiencies include the sub-decay branching fractions. Upper limits are calculated at 
the 90% confidence level and include systematics. 





Mode 


Yield 


Efficiency (%) S (xl0- 




snf. 


B+ 
50 - 
B - 




or: c+lO.l +5.6 
•Jd.d _g 4 _3 9 

+5,1+1.6 
' -4.3 -1.5 
.t; A +11.2 +7.4 
^O.U _;^o.4 -5.1 


3.40±0.20 
0.88 ±0.12 
4.28 ±0.25 


f5 q+2.0+1.2 
^••^ -1.8 -0.9 

7 Q +3.8 +1.6 
' -"J -3.2 -1.5 
f; n +1-7 +1.3 
-1.6 -1.0 


4.3 
2.5 
5.0 


B+ 


i^3*(1780)+7 


< 21.5 


0.77 ±0.34 


< 33 






50 - 


K3* (1780)07 


< 10.0 


0.17 ±0.08 


< 72 






B - 


.^1(1780)7 


< 27.3 


0.94 ±0.41 


< 34 
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